Electron localization by a magnetic vortex 
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Abstract 



We study the problem of an electron in two dimensions in the presence 
of a magnetic vortex with a step-like profile. Dependending on the values of 
the effective mass and gyromagnetic factor of the electron, it may be trapped 
by the vortex. The bound state spectrum is obtained numerically, and some 
limiting cases are treated analytically. 
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Two-dimensional electronic systems with magnetic vortices have become a topic of great 
interest, both theoretical and experimental. In the theory of the fractional quantum Hall 
effect, it is common to resort to the composite fermion (electron + magnetic flux) approach 
a similar approach has been used in anyonic systems 0. In both cases, one is left, 
therefore, with a problem of fermions moving in a random magnetic field generated by flux 
tubes. On the experimental side, random magnetic fields have been realized by using a dis- 
ordered type-II superconductor in a uniform magnetic field as the substrate for a 2D electron 
gas: the induced flux lines are trapped by the defects in the superconductor, generating a 
random array of vortices . A demagnetized ferromagnet with randomly oriented magnetic 
domains has also been used as a substrate 0. Therefore, the study of the influence of one or 
more vortices on the spectrum of a two-dimensional electron gas is a very important issue, 
especially with respect to the eventual presence of localized states. 

In this work, we study the problem of an electron in two dimensions in the presence of a 
magnetic flux (a vortex). The system is described by a Pauli Hamiltonian, and the vortex 
has the simplest profile, i.e., B{r) = By6{a — r), where 6{x) is the Heaviside step function 
and a defines the radius of the vortex. The effective gyromagnetic factor may differ from 
g = 2. In fact, although the Dirac equation predicts this value, QED radiative corrections 
give for the electron g-iactoi a value slightly bigger than 2 {qqed = 2.0023), a result 
confirmed (in fact, discovered) by experiment 0. In solids, the gf-factor of the electron 
may depart significantly from 2 as a result of spin-orbit coupling effects (for instance, 
g = 0.32 in n-GaAs and —44 in n-InSb). However, as we shall see shortly, what determines 
the existence of bound states is the quantity gm* /m = g*, where m* is the effective mass of 
the electron {g* = 0.026 in n-GaAs and —0.66 in n-InSb 0) 

The Hamiltonian has the following form: 

where fis = eh/2mc is the Bohr magneton, and cr^ = diag{l, —1) is the third Pauli matrix. 
The magnetic field points in the z-direction. In units such that 2m* = h = c = 1, this 
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simplifies to 



H = (-zV - eA)' - ^g*eBa,. (2) 



In the case of a magnetic field with radial symmetry, one can choose the vector po- 
tential in the symmetric gauge, A = A{r) e^, so that the magnetic field is written as 
5 = i A [ry4(r)], and Eq. (H) becomes 

H = --—(r—]-- — + 2ieA - — + e^A^ - - g*eBa, (3) 
r dr \ dr J dcp'^ r d(j) 2 ^ 

Thus, writing the eigenf unctions of H as '?/'(r, cj)) = Ui{r) e*^"^, we obtain the eigenvalue 
equation 

f-T^--T- + 4- 2^^^ - + - I 9*eBa)i u,(r) = E w,(r), (4) 

y dr^ r dr r'^ r 2 J 

which is equivalent to the pair of equations 

The — (+) sign corresponds to spin-up(down) electrons. 

Now, let us consider the field of a vortex, B = BvO{a — r). The vector potential 
corresponding to this field is 

Byr for r < a, 

(6) 

Bya /2r for r > a. 

Since we are interested in the possible existence of bound states, we shall consider only the 
spin-up sector of Eq. (in the spin-down sector, the coupling of the spin with the field 
gives rise to a repulsive potential). For E_ = — /i^ < 0, this equation has the following 
solution (inside and outside the vortex of radius a) 

uUr) = Q„e-*^'/2a^ {r/af\ M {a, \i\ + 1, ^rVa^) , (7a) 
Uout{r) = CoutK\i^^\{^r/a), (7b) 

where $ = eByo^ /2 is the magnetic flux in units of the quantum flux (here we assume $ 
positive), ^ = fia, 



4$ 2 V 2 
M{a, b, z) is the Kummer (or confluent hypergeometric) function and Ku{z) is the modified 
Bessel function 0. 

In order to find the energies of the bound states, one must impose continuity of the 
logarithmic derivative of u{r) at the boundary of the vortex, or, equivalently, look for the 
zeroes of the function 



Uoutir) Uinir) 



^1 ^-^W(0 ^ M(a + l,K| + 2,$) 

Now, we shall analyze (Q) numerically for representative values of the parameters, and 
consider the analytical solution of some limiting cases to provide some physical insight. 

For (yf* = 2, a theorem due to Aharonov and Casher |T0| tells us that there exists zero- 
energy normalizable states, where is an integer defined through $ = A^ + e, 0<e<l. 
Let us show this explicitly: taking ^ — > in Eq. (^, the last term vanishes if £ > 0, and the 
second term tends to — 2 |£ — $|. Therefore, 

i?(0,$,^,2) = -|^- $1 + $ (10) 

which vanishes for £ < $. There is a further restriction on the values of C.. When ^ = 0, the 
solution outside the vortex, Eq. (|7^), must be replaced by 

Uout{r-i = Q) = Coutr-\'-''^. (11) 

In order for ui{r) to be square-integrable, we need |^ — $| > 1. Together with the condition 
< £ < $ and the definition of A^, we obtain the following inequality: 0<£<$ — 1<A^. 
Thus, except for a normalization factor, the corresponding radial wave functions are 

, , f e"*''^/^"^ (r/aV for r < a, 

u^f\r) = \ (£ = 0,l,...,Ar-l) (12) 

[e */^(r/a)^ * for r > a. 

A remark is in order here: such states can only exist if a 7^ 0. For a ^ 0, they develop a 
non-integrable singularity at the origin. 



If g* > 2, we can find eigenstates with negative energy. In order to create a bound state 
with angular momentum i, a minimum value of $ is required; let us call it $^ (clearly, for 
$ > this bound state will be pushed to negative energy). To find its dependence on i, 
one must solve the equation -R(0, i, g) = 0, which, for £ > 0, can be written as 



0, 



(13) 



+ 1 M(ao, £+!,$£) 

with ao = = 0) = (2 - g*)/A. For £ = 0, it has $o = as solution. For i > 0, Eq. (|13D 
is too complicated to be solved analytically, but it becomes tractable if (7* — 2 ^ 1. In 
fact, assuming that $£ = £ — 5, with 5 <^ 1, and expanding Eq. (0) to first order in small 
quantities, one finds 

g*-2 



M(l,£ + 2, 



(14) 



4(£+l) 

Note that there is a discontinuity at g* = 2, which is not apparent in the expression above. 
In fact ^i{g* = 2) = and not i, as sugested by the g* ^ 2 limit of (p!^). This happens 

because even for an infinitesimal binding energy the wavefunction has an exponential decay 
outside the vortex, thus eliminating the requirement that |£ — $| be greater than one to 
guarantee the normalizability of the wavefunction. The behavior of as a function of g*, 
for different values of angular momentum is shown in Fig. 1. 

Now, let us examine how the energy of the bound state with i = depends on the 
magnetic flux. Assuming that ^ <^ 1 for $ ^ 1, we can approximate R{^, $,£ = 0,g*) by 



2$ <^ 1 



r(i-$) f^' 



2<S>' 



- "0 



r(i + <i>) \2j 

In this approximation, solving R{^, $, 0, g*) = for ^ we find 



l--(l-ao)$ 



(15) 



^{<l>,i = 0,g*)=2 



r(l-$) [g* + 2 [l + l(g*-2)^^ 



1/2* 



(16) 



Since we are assuming that $ ^ 1, we can write r(l ± $) ^ 1 =1= 7$ (where 7 = 0.577 ... is 
Euler's constant) and use lima;^o(l + x)^^^ = e to simplify (p!6|) to 



e($,0,,*)-2e--^V4 [^^^ 



l/2'I> 



(17) 



The full spectrum may be obtained numerically for particular values of g*, for different 
values of the angular momentum. Illustrative results are presented in Fig. 2. For the sake of 
comparison, the approximate result for ^($,0,(7*), Eq. ([T^), is plotted with the exact result 
in Fig. 3. 

For g* > 6, Eq. (^) may posses more than one solution, each one corresponding to a 
different bound state. This happens because ao + 1 becomes negative, so that M(ao + 1, ^ + 
2, has at least one non-zero root |Tl|; if such a root is larger than i, it is also a solution 
of Eq. (plsl) . Since, as far as we know, such large values of g* are not found in Nature, we 
shall not analyze them further. 
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Figure Captions: 

Figure 1: vs. g* for £ = 1 (solid line), ^ = 2 (dashed line) and £ = 3 (crosses). (For 
g* = 2, see remark in the text.) 

Figure 2: ^ vs. $ for £ = (solid line), i = 1 (dashed line), and i = 2 (crosses): (a) g = 2.1, 
(b) g = 2.0023. 

Figure 3: ^ vs. $ for g = 2.1 and £ = 0: exact (solid line) and approximate, as given by 
Eq. (|13) (dashed line). 
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Figure 1 




0.5-- 



O2 2.5 3 3.5 4 4.5 5 

g 



1- 
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